Phase diagram of depleted Heisenberg model for CaV409 
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We liave numerically investigated the 1/5-depleted Heisen- 
berg square lattice representing CaV409 using the Quantum 
Monte Carlo loop algorithm. We have determined the phase 
diagram of the model as a function of the ratio of the two dif- 
ferent couplings: bonds within a plaquette and dimer bonds 
between plaquettes. By calculating both the spin gap and the 
staggered magnetization we determine the range of stability 
of the long range ordered (LRO) phase. At isotropic cou- 
pling LRO survives the depletion. But the close vicinity of 
the isotropic point to the spin gap phase leads us to the con- 
clusion that already a small frustrating next nearest neighbor 
interaction can drive the system into the quantum disordered 
phase and thus explain the spin gap behavior of CaV409. 
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The stabihty of the long range ordered (LRO) ground 
state of the planar Heisenberg model has been at the fo- 
cus of investigations for a long time. The recent discovery 
of a spin gap in CaV409 has given special importance 
to this question. This compound can be described by a 
1/5-depleted planar antiferromagnetic Heisenberg model 
. One of the important questions regarding this ma- 
terial is whether the depletion of the square lattice can 
account for the spin gap, or if additionally frustration 
effects are important. 

The role of lattice defects and depletion in destabiliz- 
ing LRO has been studied in a variety of contexts. One 
example are spin ladders which can be obtained from 
the planar copper oxide materials, by breaking up the 
planes into ladders of constant width [Q. Another way 
to destroy LRO is to deplete the lattice. The bonds be- 
tween spins are then weakened, similar to the introduc- 
tion of holes, and quantum fluctuations are enhanced, 
which might destroy LRO. An example is the triangu- 
lar Heisenberg antiferromagnet, which exhibits LRO. De- 
peletion of 1/4 of the spins, leads to the Kagome-lattice, 
which is believed to have no LRO 

The stability of LRO is also of relevance in the field of 
high temperature superconductors. There the rapid de- 
struction of LRO upon hole doping and the possibility of 
realizing a doped resonating valence bond (RVB) phase 
exhibiting a finite gap in the spin excitation spec- 
trum (spin gap) are of great current interest. The study 
of lattice defects, such as depletion or the formation of 
ladders can give valuable insights 0. 

The lattice structure of CaV409 and the 8-spin unit 
cell [|| used in our simulations is shown in Fig. |l|. It 
can be viewed as consisting of loosely connected 4-spin 



plaquettes. Two topologically different types of bonds 
can be distinguished. One are bonds within a plaque- 
tte Jo, the other are dimer-bonds connecting plaquettes 
Ji. Additionally CaV409 is believed to have a significant 
frustrating next nearest neighbor (n.n.n.) antiferromag- 
netic interaction 

Ueda et al. |Q| and Imada and Katoh |^ have argued 
that the spin gap can be explained as originating in a 
plaquette RVB state, consisting of local singlets of the 
four spins on a plaquette. This plaquette RVB state is 
the exact ground state in the limit Ji = 0. Second or- 
der perturbation theory around this limit suggests that 
it survives even at isotropic coupling Pj^]. A perturba- 
tion around the dimer limit Jo = [0 also leads to a 
wide range of stability of the dimer singlet phase, but 
a small range of the couplings exists, where no gap is 
observed in perturbation theory. First Quantum Monte 
Carlo (QMC) results by Katoh and Imada ||] also suggest 
the existence of a finite gap A = 0.11±0.03at isotropic 
coupling. 

Linear spin wave theory (LSW) § and Schwin ger bo- 
son mean field theory results jl^] on the other hand indi- 
cate that LRO could survive at isotropic coupling despite 
the depiction of the lattice. Exact diagonalization results 
are also contradictory 1 1 1 . They suffer greatly from 
the restriction to very small clusters and the extrapo- 
lation to the infinite system size is difficult. Sano and 



Takano 1 1 1 1 and Albrecht and Mila [|10[ find a small spin 
gap, but also a substantial staggered magnetization [ p^ . 
No definite conclusions can thus be drawn from these 
calculations either. 

To resolve these conflicting results we have determined 
the phase diagram (see Fig. ^ of the non-frustrated 
model using the QMC loop algorithm ||l2|. Using this 
highly efficient cluster method we can investigate larger 
systems at lower temperatures and with a much higher 
accuracy than possible with the standard world line al- 
gorithm used by Katoh and Imada ||^. We have investi- 
gated lattices with up to = 800 spins at temperatures 
down to T = 0.02. The QMC method suffers from no 
systematic errors and the results are reliable within the 
statistical errors. 

We find a wide region of stability of the Neel-ordered 
phase as a function of the ratio of the couplings a — 
Jq/Ji- We estimate the lower boundary to lie between 
0.55 < < 0.65 and the upper boundary between 
1.05 < < 1.1. At isotropic coupling LRO thus sur- 
vives the depletion of the lattice. The critical point a" is 
quite close to isotropic coupling, and a small frustration 
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might be sufRcient to drive the system into the disordered 
state. 

To determine this phase diagram we have calculated 
both the spin gap and the staggered magnetization. The 
spin gap A can be obtained from the low temperature 
behavior of the uniform susceptibility x- Figure ^ shows 
x(r) for some representative points. In a gapped system 
it decreases exponentially as e~^/^ for low temperatures. 

Any finite system exhibits a gap, and thus a careful 
treatment of finite size effects is necessary. For each tem- 
perature we have done calculations on clusters of differ- 
ent size (up to = 800 spins) to see whether our results 
have converged to the infinite system size limit. In the 
regions of a large gap the convergence is quite rapid and 
it is no problem to obtain the gap A from a fit of the low 
temperature behavior of the uniform susceptibility x to 
an exponential decay e~^/'^ . In case of a vanishing or 
very small gap on the other hand the susceptibility de- 
creases linearly down to the lowest temperatures Tq we 
could study rehably on our finite clusters. In these cases 
we cannot definitely decide about the existence of a gap, 
but can only give an upper bound A < Tq. 

In Fig. H we plot the gap obtained in this way together 
with the perturbation theory estimates [§,1). Perturba- 
tion theory is surprisingly accurate, but overestimates 
the gap slightly. Specifically at the isotropic point we 
do not see any indication for a gap, in contradiction to 
Katoh and Imada Their calculation of x(T) is for 
a much smaller lattice (80 spins), and their gap may be 
due to finite size effects p^ . 

The existence of LRO can be checked by calculating 
the staggered magnetization m^: 
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•nis vanishes in the infinite system size limit in case of 
purely short range correlations, while it is finite for LRO. 
The finite size scaling of is known and a reliable ex- 
trapolation possible uM : 



ms{N) = TOs(oo) + O 



1 



(2) 



Figure ^ shows the system size dependence of m^. Let 
us first discuss couplings in the spin gap regime. There 
the finite cluster results can be extrapolated linearly in 
1 / \/N to zero moment in the infinite system [Fig. ^(a)] . 
In the double logarithmic plot [Fig. ^b)] it can clearly 
be seen that the results for finite clusters bend down and 
approach the linear decrease (slope 1). The results for 
couplings in the LRO phase on the other hand bend up 
and reach a constant value asymptotically. At the criti- 
cal coupling itself we expect a power law with a critical 
exponent different from the ^^-behavior of the gapped 
phase. A rough estimate shows an exponent of the or- 
der 0.5, as expected from the mapping to the non-linear 



a model [|l5[. This exponent is indicated as a dotted 
line. A more detailed investigation, to obtain a reliable 
estimate for the exponent and a better estimate for the 
critical coupling is currently under progress. 

Although the system size dependence is asymptotically 
linear in N^^^"^, our lattices are not yet large enough to 
be really in that limit. To get an estimate for the qual- 
ity of our extrapolations we extrapolate both rus and 
m^. In case of LRO both extrapolations should be lin- 
ear. We observe that, as seen in Fig. |^(a) the system 
size dependence is not perfectly linear, but still bends 
down a little bit. Thus we take the value obtained from 
this fit as an upper bound. In a plot of mj. on the other 
hand a slight upwards bend can be observed and we take 
that extrapolation as a lower bound. Both extrapola- 
tions agree well. In the phase diagram (Fig. ||) we show 
the average value, with the error bars indicating these 
upper and lower bounds. Wc have tested this procedure 
for the square lattice, where our result of nis — 0.306(3) 
agrees perfectly with the most accurately known value 
TO, = 0.3074(4) |l|]. A gain close to the critical points 
the moment is very small and a definite decision about 
a nonzero magnetization difficult. The magnitude of the 
staggered moment compares well with the results of lin- 
ear spin wave theory (LSW) (also shown in Fig. but 
the range of stability of the LRO phase is overestimated 
by the LSW. 

The conclusions obtained from the estimation of the 
gap and the staggered moment are perfectly consistent. 
Starting from the dimer limit we see a decrease of the 
gap as Jq is increased. At a = Jo/Ji = 0.55 we can still 
find a finite gap, while at a = 0.65 we observe a finite 
staggered magnetization and a zero or small gap. Thus 
we conclude that at a critical coupling 0.55 < a[. < 0.65 
the dimer singlet phase becomes unstable and the model 
exhibits LRO. The critical coupling is probably close to 
a = 0.6. There we cannot definitely decide about the 
existence of a gap or LRO from our finite cluster results. 
Starting from the plaquette side the gap also decreases 
as we increase Ji, but the plaquette RVB state is sta- 
ble for a wider range of couplings than the dimer state. 
This is quite natural, as each spin is connected to one 
dimer bond, but two plaquette bonds. Perturbation the- 
ory predicts that the isotropic point is still in the range 
of stability of the plaquette RVB state, but our QMC 
simulations show that LRO sets in at 1.05 < a" < 1.11. 
At the isotropic point we observe a substantial nonzero 
staggered magnetization — 0.178(8). 

Comparing our results to previous calculations we find 
that the region of stability of LRO is larger than es- 
timated by second order perturbation theory [^, but 
smaller than estimated by linear spin wave theory and 
Schwinger boson mean field theory [|lO| . Our results also 
agree well with the exact diagonalization estimates of the 
staggered magnetization jlOj, while the extrapolation of 
the spin gap data by exact diagonalization is unreliable 
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By varying the ratio of the couphngs in the 1/5-th de- 
pleted square lattice we can study both the LRO phase 
and the disordered phase of a two dimensional quantum 
antiferromagnet, without having to introduce frustration 
or to break symmetries, as in the dimerized square lat- 
tice model iQ. This model is thus ideal to study the 
critical behavior and to test the predictions made by 
Chakravarty, Halperin and Nelson based on the 2 + 1- 
dimensional non- linear a model [ p^ . 

In comparison to experimental results on CaV409 we 
conclude that the depletion of the square lattice alone 
is not sufficient to destroy LRO in the Heisenberg an- 
tiferromagnet, but it is very close to the critical point. 
An additional frustrating next nearest neighbor coupling 
is needed to drive the system into the gapped plaquette 
RVB phase. All estimates from perturbation theory ||] 
and exact diagonalization |ll] agree that the stability of 
this plaquette RVB phase and the gap are greatly en- 
hanced by a frustrating next nearest neighbor coupling. 
Thus we expect that already a quite small frustration 
will destroy LRO and can explain the substantial gap 
observed in CaV409. 
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larger than the nearest neighbor coupling. This question is 
of no relevance to the validity of our results here, since the 
lattice obtained by considering only the n.n.n. coupling are 
just two non-interacting copies of the same lattice. The in- 
clusion of the n.n. coupling on the other hand would then 
lead to a different frustrating coupling, but the main effect, 
frustration would be the same. 
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FIG. 1. The lattice structure of the depleted Heisenberg 
lattice describing CaV40g. Indicated are the two different 
types of bonds, plaquette-bonds Jo and dimer bonds J\. 
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The smallest possible unit cell contains 4 spins. For tech- 
nical reasons we have chosen an 8-spin unit cell in our 
simulations. 

There are still some discussions going on about the mag- 
nitude of the couplings. It has been argued (see e.g. Ref. 
[|ll| that the next nearest neighbor coupling could even be 



FIG. 2. Phase diagram as a function of the ratio Jo/ Ji. (a) 
shows the whole range of couplings. The leftmost point cor- 
responds to the dimer limit Jo — and the rightmost point 
to the plaquette limit Ji = 0. (b) A detail of the phase dia- 
gram around the isotropic point plotted as function of Jo/Ji. 
Circles indicate our QMC results for the spin gap, normal- 
ized by Jo + Ji. In the gap less region the error bar indicates 
an upper limit for the gap. Diamonds show the staggered 
magnetization. The error bars indicate the upper and lower 
bound, as described in more detail in the text. As reference 
we have included the perturbation theory estimates for the 
gap Q and the linear spin wave theory (LSW) estimates for 
the staggered moment. 



FIG. 3. Temperature dependence of the uniform suscepti- 
bility X for different ratios of the couplings Jo/Ji. For each 
temperature the system size was taken large enough to see the 
value for the infinite system. The lowest temperatures were 
calculated on a A^ = 512 spin lattice. As a reference we have 
included results for the square lattice Heisenberg model. The 
temperature is in units of the larger of the couplings Jo and 
Ji. 
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FIG. 4. System size dependence of the staggered magne- 
tization rus for different ratios of the couphngs Jo/Ji- For 

each system size the temperature was chosen low enough to 
see the ground state properties. The largest systems contain 
N = 800 spins, (a) ma plotted as a function of N~^^^. A 
linear extrapolation gives the bulk value, (b) A double log- 
arithmic plot clearly shows the existence of long range order 
or the linear decrease with system size respectively. Included 
as guides to the eye are two straight lines corresponding to 
power law decays with powers 1 and 0.5. 
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Figure 2, Troyer et al. 
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Figure 4, Troyer et al. 





